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Using the complex Toda chain (CTC) as a model for the propagation of the iV-soliton pulse 
trains of the nonlinear Schrodinger (NLS) equation, we can predict the stability and the asymptotic 
behavior of these trains. We show that the following asymptotic regimes are stable: (i) asymp- 
totically free propagation of all N solitons; (ii) bound state regime where the N solitons move 
quasi-equidistantly; and (iii) various different combinations of (i) and (ii). We show with exam- 
ples of N = 2 and 3 how this analysis can be used to determine analytically the set of initial 
soliton parameters corresponding to each of these regimes. We also compare these analytical re- 
sults against the corresponding numerical solutions of the NLS and find excellent agreement for 
each of the regimes described above. We pay special attention to the regime (ii) because the 
quasi-equidistant propagation of all N solitons is of importance for optical fiber soliton commu- 
nication. Our studies show that such propagation can be realized for N = 2 to 8. For the case 
of 3-soliton trains we also show how one can determine the instability regions and the transition 
regions between the various regimes. Finally we propose realistic configurations for the sets of 
the amplitudes, for which the trains show quasi-equidistant behavior to very large run lengths. 
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I. INTRODUCTION 

One of the important problems in optical fiber soliton 
communication is to achieve as high of a bit rate as pos- 
sible. In order to do this, one needs to be able to pack 
the solitons into as short of a space as possible. However, 
if the solitons are too close together, then their mutual 
(linear and nonlinear) interactions can cause them to col- 
lide and/or separate, thereby corrupting the signal. The 
current solution of this problem is simply to require each 
soliton to be sufficiently far apart from all others (usu- 
ally 6 or so soliton widths) so that such interactions can 
be totally neglected. However, at the same time, it was 
predicted and experimentally confirmed |s| that for 
certain values of relative soliton parameters, this separa- 
tion can be reduced, and at the same time, still maintain 
signal integrity. Our main purpose here is to analytically 
and numerically detail the soliton parameter regime, in- 
side of which, signal integrity can be maintained. In par- 
ticular, we are interested to determine how one may use 
this inter-soliton interaction for stabilizing a soliton train. 

Any optical communication signal will be composed 
of "random" combinations of O's and l's. Such a sig- 
nal can also be viewed as being composed of a random 
collection of iV-soliton trains, with varying widths of O's 
between them. Thus it is then adequate for us to simply 
analyze the stability of individual iV-soliton trains, for 
iV = l,2,3,.... This we will do and will study a nonpe- 
riodic and finite train (chain) of soliton pulses by both 
analytical and numerical methods using and developing 
the ideas in [j||J. 

The basic model and description of 7V-soliton trains in 
optical fibers is provided by the nonlinear Schrodingcr 
(NLS) equation and its perturbed version: 



iu t 
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|zi| 2 u(x, t) = iR[u]. 



(i) 



This equation describes a variety of wave interactions, 
including solitons in nonlinear fiber optics [l|J^-|l5t| and 
spatial solitons in nonlinear refractive media pi[ ■ 

The inverse scattering method HQ allows one to solve 
exactly Eq. (Q) when R[u] = and to calculate explicitly 
its A^-soliton solutions. However for our purposes, this 
method is impractical for two reasons. First, there will be 
important physical problems in this system in need of ad- 
dressing where R[u] will be nonzero, for which no explicit 
solutions are known. Second, an approximate method 
can serve much better than an exact approach since the 
7V-soliton trains that we need to study are rather special 
and can only be approximated by A^-soliton solutions. 
Such trains are actually sums of I-soliton pulses, which 
are spaced almost equally, have almost equal amplitudes, 
and move with essentially the same velocity. More specif- 
ically, they are the solutions to Eq. ([!]) with R[u] = 
satisfying the following initial conditions: 



u(x, 0) 



N 
fc=l 



«(a:,0), 



(2) 



where Ui a -k(x, t) is the 1-soliton solution of the NLS given 
by: 



U la ;k(x,t) 



2v k e i<1 "> 



cosh(2^ fe (x - &(t)) ' 



8 k {t)=2{y 2 +v 2 )t + 5 kfl 



(3a) 



(3b) 



(3c) 



(3d) 



Here by S k (t), £&(£) we have denoted the phase and posi- 
tion respectively of the fc-th soliton, 2v k is soliton's am- 
plitude while is its velocity. The quantities S k $, £ k Q, 
i^fco, and /ifco all denote the initial values at t = 0. 

An effective method for studying the interaction of 
such trains of soliton pulses was first proposed by Karp- 
man and Solov'ev (KS), for the simplest non-trivial case 
of a 2-soliton interaction |l6| . Further developments and 
analysis for different physically important perturbations 
can be found in Jl7],|l8| and the references therein. For 
alternative approaches see e.g. |f9| p2| : for a review see 
Ref. §. 

The KS method is based on the adiabatic approxima- 
tion. It is valid for any collection of well separated soli- 
tons, such that their mutual interactions will lead to a 
slow deformation in the soliton parameters. These con- 
ditions are met provided the soliton parameters satisfy: 



|MfeO — MnOI "C MOj 
K-0 - V n o\ < VQ: 
^olCfeO - £no| > I, 
K-0 - ^no||6cO - 61 I < 1, 



(4a) 
(4b) 
(4c) 
(4d) 



where 2v$ and 2/io are the average initial amplitude and 
velocity. 

How fast adjacent solitons will interact can be mea- 
sured by a control parameter, eo, which is a measure of 
both the overlap between the neighboring solitons and 
also the strength of their interaction. This is given by 
eo = h , oe~ 2l/ ° r ° , where ro is the average initial pulse sep- 
aration. 

Gorshkov and Arnold |24j ], have conjectured that 
an infinite train of out-of-phase soliton pulses, with equal 
amplitudes and velocities could be described by the real 
Toda chain: 



d 2 q k 
dr 2 



olk+l—Qk 



(5) 
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where r = ^v^t, k — 0, ±1, ±2, . . . and qk are real func- 
tions related to the soliton positions. This system will 
be referred to as the real Toda chain (RTC). The next 
step towards the physically more realistic case with fi- 
nite number of solitons was proposed in p5[ . In it the 
RTC was shown to describe the propagation of N equal 
amplitudes out-of-phase solitons numerically. 

Recently in Refs. |JJ^,^,^|, the Karpman-Solov'ev 
method was extended to A-soliton pulses, and then with 
additional approximations, was reduced to the complex 
Toda chain equations (CTC) (||) with N sites. The cor- 
responding system of equations is (||), but with k — 
l,...,N and e~ qo = e qN+1 = 0. Meanwhile, the complex 
valued functions qk(t) are related to the soliton parame- 
ters by: 



qt+i — qk — — 2fo(6c+i — £fc) 

+ i (tt + 2/x (£ fc+ i - £fc) 



In Avl 



(S k+1 - S k )) . (6) 



With this, the problem of determining the evolution of 
an NLS A-soliton train has been reduced to the prob- 
lem of solving the CTC for TV sites. Since (||) is also 
integrable, then we may use the special techniques valid 
for integrable lattices (and chains) to study this prob- 
lem. It has already been shown that one may determine 
the asymptotic behavior of asymptotically separating A- 
soliton trains, simply by analyzing the eigenvalues of a 
certain matrix 0]. 

Our main results, briefly reported in |27| are the fol- 
lowing: 

(i) we show by analyzing the exact analytic solutions 
of CTC, that it has several qualitatively different classes 
of asymptotic regimes. Besides the asymptotically free 
motion (which is the only possibility for the RTC) , CTC 
allows also for: (a) bound state regime when all the A 
particles move quasi-equidistantly; (b) all possible inter- 
mediate regimes when one (or several) group(s) of parti- 
cles form bound state(s) and the rest of them go into free 
motion asymptotics. In addition to these relatively stable 
regimes of motion there are also less stable regions in the 
space of soliton parameters, where one regime switches 
into another one. There one can find (c) singular solu- 
tions, which tend to infinity for finite values of r, and 
(d) various types of degenerate solutions when two or 
more of the eigenvalues become equal. 

(ii) we show, by comparing the predictions from the 
CTC model with the numerical solutions of the NLS, 
that regimes of type (a) and (b) indeed take place in the 
soliton interactions and are very well described by the 
CTC model. Our analytic approach allows us to predict 
the set of initial parameters, for which each of the asymp- 
totic regimes takes place. We put special stress on the 
bound state and quasi-equidistant regimes (a) since such 
solutions would be optimum in long distance fiber optics 
communications. 



II. ASYMPTOTIC REGIMES OF THE CTC 

As in one can generalize the RTC p9|-p3[| to the 
complex CTC case. We list the four most important 
points concerning this below: 

SI) The CTC Lax representation is the same as for the 
RTC: 



L=[B,L], 



N 



L — ^ {bkEkk + ak{E k ,k+i + : 



fe=l 



a? 



B — ^ dfc (EkM+i — Ek+i, k ) 



(7a) 



(7b) 



(7c) 



k=l 
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Here the matrices (Ek n ) p 
ever one of the indices becomes or A 
notations in (J?!) are as follows: 



,(«k+i-«*)/ a 



when- 
1: the other 



(8a) 



b k = - (jUfc + iVk) 



(8b) 



52) The integrals of motion in involution are provided 
by the eigenvalues, C,k, of L. 

53) The solutions of both the CTC and the RTC 
are determined by the scattering data for Lq = L(r = 
0). When the spectrum of Lq is nondegenerate, i.e. 
Cfe 7^ Cj f° r k j, then this scattering data consists 
of {Cfc) r fc}fcLi- Here are the first components of the 
corresponding eigenvectors ^ fc - > of Lq: 



Lot 



(*0 - 



(fc) 



(9) 



which are uniquely determined (up to an overall sign) by 
the normalization condition 



N 

E(e. 



(fc) 



(10) 



see 



(30 -32]. This scattering data uniquely determines 
both Lq and the solution of the CTC. 

S4) Lastly, the eigenvalues of Lq uniquely determine 
the asymptotic behavior of the solutions of the CTC; 
these eigenvalues can be calculated directly from the ini- 
tial conditions. We will extensively use this fact for the 
description of the different classes of asymptotic behav- 
ior. 

In addition to the dynamical variables becoming com- 
plex valued, there are other, important differences be- 
tween the RTC and CTC, and their spectra. For the 



RTC, one has that |30|j3 1|] both the eigenvalues, Ck, and 



the coefficients, r^, are always real- valued. Moreover, one 
can prove that Cfc 7^ Cj f° r k 7^ j, i.e. no two eigenvalues 
can be exactly the same. As a direct consequence of this, 
it follows that the only possible asymptotic behavior in 
the RTC is an asymptotically separating, free motion of 
the solitons. 

This situation is different for the CTC. In addition 
to all dynamical variables being possibly complex, the 
eigenvalues can also become complex, = Kk + *?7fc, as 
well as the coefficients, rfc, as given above. Furthermore, 
one could now have multiple eigenvalues. However, the 
collection of eigenvalues, still determines the asymp- 
totic behavior of the solitons. In particular, it is Kk that 
determines the asymptotic velocity of any soliton. For 
simplicity, here we shall always take Cfc 7^ Cj f° r k ^ j. 
(However, this condition does not necessarily mean that 
Kfc 7^ Kj.) We also may assume that the k^s are or- 
dered as: K\ < «2 < . . . < K.N- Once this is done, then 
in any train of solitons, there are three possible general 
configurations: 

Dl) Kk 7^ Kj for k 7^ j. Since the asymptotic velocities 
are all different, one has the well known asymptotically 
separating, free solitons. 

D2) Ki — K2 — ■ ■ ■ — kn. In this case, all N solitons 
will move with the same mean asymptotic velocity, and 
therefore will form a "bound state". The key question 
now will be the nature of the internal motions in such 
a bound state. In particular, one would want any two 
adjacent solitons to move quasi-equidistantly. 

D3) One may have also a variety of intermediate situa- 
tions when only one group (or several groups) of particles 
move with the same mean asymptotic velocity; then they 
would form one (or several) bound state (s) and the rest 
of the particles will have free asymptotic motion. 

Obviously the cases D2) and D3) have no analogies in 
the RTC and are qualitatively different from Dl). The 
same is also true for the special degenerate cases, where 
two or more of the £fc's may become equal. These cases 
will be considered elsewhere. 

Another type of solutions of the CTC which should 
be dealt with separately are the singular solutions. They 
appear in a special non-euclidean formulation of the RTC 
J33|; they can be obtained from CTC by taking one or 
several of the a& 's purely imaginary and the others - real- 
valued functions of r. 

In order to avoid complicated and long formulas, we 
will skip most of the technical details and will limit our- 
selves with the simplest nontrivial cases: N = 2 and 
N = 3. We also assume without loss of generality that 
tr L = T,k=i Cfc = °- 



A. The N = 2 case 

The general solution of N — 2 CTC can always be 
chosen so that qi(r) = — ^(t) (the center-of-mass is at 



rest and at the origin of the coordinate system). Then 
(see also [pi): 



<?i(t) = In 



cosh(2Cir + 7) 
20 



where 7 = — ln(ri/ra). From formula (Q) we get: 

/x , s , „ 7r — S(t) 
qi(r) = v Q r( T )-\n2v + i -A-L, 



(11) 



(12) 



where r(r) = &(r) — £i(t) and S(t) — S2 (r) — 5i(t) 
are respectively the distance between the solitons and 
their phase difference. From @, (0) and (||), we find the 
following expression for L$, in terms of the initial soliton 
parameters: 



-3 (A/Lio + iAv ) iv Q e- Uor °- i&a / 2 
iv e- Voro - lS °/ 2 i(A/i,^A^) 



(13) 



where A^ = ^ 20 - fiio, Av = is 2 o - vvs and S =5(0). 

For simplicity, from now on, we will consider only those 
trains where the solitons have vanishing initial velocities: 
/ifco = in some moving coordinate system. Then solving 
the characteristic equation for Lq, we find: 



Ci,2 = = ±«i ± ir/l, 



D = (Az/ ) 2 + 16i$e 



-2v r -iSo 



(14a) 



(14b) 



When D is complex in general, with a nonzero imagi- 
nary part, then by the above, n\ will be nonzero and the 
two solitons will always asymptotically separate. How- 
ever, when D is real, then it is possible for a bound state 
of two solitons to form. This requires the relative phase 
between the solitons, So, to be either or w. If it is zero, 
then K% is always zero, regardless of the size of any am- 
plitude variations. If it is n, then k± will again be zero, 
provided that \Ai/q\ > Az/ cr! 2, where 



Au 



cr,2 



4i/ e 



(15) 



and 2Av cl _2 is a critical value for amplitude variations in 
a 2-soliton train. This means that if two adjacent solitons 
have a sufficiently large initial difference in their ampli- 
tudes of 2A^o > 2Ai/ crj 2, and if the phase difference is 7r, 
then these two solitons will never asymptotically sepa- 
rate. In other words, a slight variation in the amplitudes 
can stabilize and prevent two solitons from asymptoti- 
cally separating. Obviously, in this case the distance, 
r(t) (r = Avot), between the solitons will be periodic in 
t with period: 



1 2 — 



2^ 0V /(A^) 2 ±(A^ 2 ) 2 



(16) 



where the signs plus and minus correspond to Sq = and 
So = 7r respectively. 
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Most importantly, this motion is then bounded, and 
optimally its range should be significantly smaller than 
the initial spacing between the solitons. 

Define 



A k = 



^fc,max T fe,min 

ro 



(17) 



where f fcimax = max(^ + i - £ k ) and f fe , min = min(£ fe+ i - 
£k) are the maximal and minimal values of £,k+i{i~) — 
£k(r). When ^4^ <C 1, we will call such motion "quasi- 
equidistant" . For such motion, the solitons will not 
asymptotically separate, but instead will slowly oscillate 
with some small amplitude. 

For N = 2 from the explicit form of the solution, we 
can easily recover (see also [p6)): 



_ J_ z/ 2 (cosh(2 7o )-l) 
ri ' min " 2^ 2r,l 



_ 1 i/g(cosh(2 7 o) + l) 
ri, max - m ^ 



Therefore A\ equals to: 



A 1 = ^-ln C08h ^i + 1 . 
2^o?*o cosh(27o) — 1 



where 



1 y/e- u ° +y$-y 
7 = 7o + Hi = o ln / - A =4 ■ 



(18a) 



(18b) 



(18c) 



(19) 



and yo = |Ai/ |/Ai/ cr)2 . 

Let us look at our two cases. For Sq — 0, one obtains 
rimax = r and: 



average amplitude 2^o; this would diminish Av CT ^2 and 
as a consequence, also A\. 

From here on, we will need to refer frequently to vari- 
ous sets of initial conditions for the 7V-soliton trains; for 
clarity and brevity we will denote them by quadruples 
N\t \2Avq • 10 2 | — , where N is the number of pulses, r is 
the distance between neighboring pulses. In most of our 
runs, listed in the tables, the pulses are initially equidis- 
tant (i.e. £fc+i.o — £fc,o = ^o), the initial phases are chosen 
to be of the form {0, 62,0, 0, 62,0, 0> • • •} an( i the ampli- 
tudes are either equal (A^o = 0) or Afo = ^fe+1,0 — Vkfi 
and such that the average amplitude equals 2vq. The 
initial velocities, 2/ifco, in all of our runs are zero. 



B. The N = 3 case 

We assume that q\ (r) + q%{T) + qz{r) = 0, i.e. the 
center of mass is fixed at the origin. This is compatible 
with tr_L = d + ( 2 + C3 = 0. Our analysis below is 
not exhaustive since we will consider only the particular 
physically important submanifold of soliton parameters 
R, which is restricted by: a) zero initial velocities of the 
solitons; b) initially equidistant solitons £fc+i,o— £fc,o = ro. 
Condition a) means that the diagonal elements in Lq are 
purely imaginary bk — idk/i; from b) there follows that 
|ai| = 1*3-2 1 for t = 0. 

It is clear that the initial soliton parameters determine 
Lq, but the regimes of propagation are determined by 
the eigenvalues of Lq. Thus we consider various possible 
combinations of the latter. 

(i) Asymptotically free propagation. Choose 
Ki < K2 < K3, r)k - generic. Then we have: 

lim ( qi (r) + 2dr) = \ ln £ , (22a) 

T^OO 3 A12A13A23 



ro 



J- ]*-*£- 
2^o Vo + 1 



(20) 



Note that for even moderate values of ro, almost any 
small, nonzero value of Avq prevents the singularity in 
the CTC model; taking 2Ai/ appropriately large (say 
0.2 for r Q = 8, see Table IV) leads to a quasi-equidistant 
motion. 

For Sq = 7r analogously we find fi m i n = tq and: 



lim (g 2 (r) 



v2 
*12 



, , 2 P2 A{ 



lim (, 3 ( T ) + 2C3T) = |ln^Mk 

T^OO ^ Al2 



where 



(22b) 



(22c) 



ro 



1 i„. «' 



2^o Vo - 1 



(21) 



Pk 



rir 2 r 3 



(22d) 



Note that in this case, the quasi-equidistant regime holds 
only for yo > 1. The value yo — 1 is a critical point, at 
which the quasi-equidistant regime switches over into the 
asymptotically separating regime. 

Thus the motion will be quasi-equidistant if e 2 ' 70 ^> 
1. Both formulas (p0[), ( |2l| ) show that an increase in 
|Az/o| diminishes the oscillations of r(r). Another way 
to diminish A\ for fixed Az/q and ro is to increase the 



A 4fe = 2(0 - CO 



(22e) 



Obviously case (i) corresponds to asymptotically free mo- 
tion of the pulses; their velocities are determined by Kk- 
This case was investigated in j|,|5| for 7V|r |0|l type 
trains. We note that for this type of initial conditions 
the CTC reduces to the RTC, for which the asymptoti- 
cally free regime is the only possible one. 
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(ii) Three-soliton bound state. Choose k% = k 2 = 
K3 = and 771 = —773, r/2 = 0. Then we find the following 
periodic solutions of the CTC: 



<Zi(t) = In 



q 2 {r) = In 



g 3 (r) = In 



where £1 = ir\ x 



2cosh(2Cir 4 


-7)- 


1- P2 








2cosh(2Cir 4 


•7H 


_ _4_ 


2cosh(2Cir 4 


-7) - 




8Ci 2 






2cosh(2Cir 4 


-7)- 


P2 



-In . 

3 V 2 



) , (23a) 



~ln(?), (23b) 



1, (P2 

3 ln ly 



V(A^ ) 2 ±(Az, cr , 3 ) 2 



^cr,3 



2V2 



v e 



73 



In 2 

n 



: 73,o 4- 273,1- 



(23c) 



(23d) 



(23e) 



(23f) 



and the signs plus and minus correspond to 5q — and 
5q = 7r respectively. Obviously, if 73,0 7^ the trajectories 
of all three particles oscillate and the distance between 
them remains bounded. The corresponding period of the 
oscillations in t is equal to 



2^V(A^o) 2 ±(A^ 3 ) 2 2i/ J7f 



(23g 



Thus we conclude that this case corresponds to a 
bound state of all three particles. 

The same property is shared also by the more compli- 
cated solution for which rjk 7^ f]j 7^ take generic values 
constrained only by rji + rj 2 + 773 = 0. As we mentioned 
this solution of CTC will be periodic only if the ratio 
(Vi ~ V2) /(vi ~ V3) is a rational number. 

In the symmetrical case, i.e. for qiir) — —qz{r), 
(72 (t) = 0, we can consider A\ i2 like in ([IT]) and find that 
A\ = A%. The quasi-equidistant regime requires again 
Ax < r . 

From the explicit form of the solution we find: 



'4. 



ri, 



1 ln z/g(cosh(7 3;0 )- 1) 
2^o r?? 



1 ln ^ 2 (cosh(7 3; o) 4- 1) 
2^o Vi 



and consequently 
Ai = A 2 = 



1 



cosh 73,0 4- 1 
2v r "' cosh 73 i0 - 1 



In 



(24a) 



(24b) 



(24c) 



In terms of the soliton parameters we have 



73.0 + ?73,i = In 



-iS 



— 2(5, 



+z' 



where zq 
gets h,ma 



|A^o|/A^ crj 3. In particular, for 5q 
= rg and: 



1 



hi-P—. 

2i> z% + 1 



?"l,min = r 

For 5q = 7r analogously we get f 1 



= tq and: 



ri, 



1 

2i7 



■ In ■ 



- 1 



(25) 
one 

(26) 
(27) 



These formulas are very similar to the ones for the N = 2 
case; note however that the value of Aj/ cr 3 is different 
from the one of Av CXt 2- Note also that the symmetrical 
case is a special one and takes place only when p 2 = 2. 
If this is not the case, then the three soliton bound state 
cannot be effectively reduced to the two-soliton one. 

Like in the N — 2 case, the motion will be quasi- 
equidistant if e 2|73 <° > 1. Both formulas (§|), @ show 
that an increase in |A^o| diminishes A±. Another way 
to diminish Ai for fixed A^o and ro is to increase the 
average amplitude v§\ this would diminish Av CT .% and as 
consequence, A\ would also. These conclusions are illus- 
trated on Fig. 1 where we have plotted the four functions 
defined by: 



Ai; fcm (Az/ ) = Ai(r = 8, Ai/ , S = 8 { Q k \v = 



), (28) 



for 5q — kn, v Q m> = and k and m take values and 
1. As one can easily see, for Auq > 0.04, A\ is always 
less than 0.1. Thus with an amplitude variation of just 
more than 8%, it is possible to have a bound state with 
a rather low value of Ak ■ 

The singular behavior of Ai-Q m for A^o — > corre- 
sponds to a singular solution of the CTC; the numeric 
solution of the NLS shows that the solitons do not col- 
lide, but come rather close to each other for the values of 
t where the CTC solution develop singularities. We will 
return to this question in Section III below. 

On the other hand the singularity of ^4i ; i m for Ai^o = 
A^ cli 3 corresponds to the fact that at this critical point 
the quasi-equidistant regime switches over into the free 
motion regime. Indeed, as we will see below, for Avq < 
A^ cr ,3 we have regime (i) and the distance between the 
solitons grows infinitely while for Ai/q > Av cr ^ we get 
regime (ii) and a possible quasi-equidistant behavior. 

(iii) Mixed regime. Choose n\ < k 2 — K3, 112 7^ V3- 
Then we have: 



,(m) _ 5+m 



lim ( gi (r) +2&t) 



In- 



Pi 



3 A13A23A12 



(29a) 



lim (q 2 (r) - Cit) 



1, A 12 A 13 
- hr — Tt 

3 A 2 3 Pl 



ln(2cos(?7r + T)), 
(29b) 



lim (q 3 (r) - &r) = 1 In Al2Al3 ^ 3 - ln(2cos(7 ? r + T)), 

T-+OO Ci Pi 

(29c) 

»? = 2(772 - m), (29d) 



where 



1 = i In ; — = 1 



r 3 A 



»ri 
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(29e) 



As it is easy to see now, the second and the third parti- 
cles asymptotically form a bound state. They both move 
with the same mean velocity and the distance between 
them is bounded (as long as Ti ^ 0) and is asymptoti- 
cally a periodic function of time. 

In the particular case when Ti = 0, we get a singular 
solution. Indeed, if 7^ is such that ^Tfe+To = (k+^)ir for 
s ome integ er fc, then the right hand sides of the equations 
(291) and (29c) become infinite. 

If the CTC model predicts collisions, this always in- 
dicates that the NLS solitons are becoming dangerously 
close together, see e.g. [p|. There is an excellent match 
between CTC and NLS except in the vicinities of the 
points where the CTC develops singularities since here, 
(4) is becoming violated. 

Let us now describe the particular choices of the soliton 
parameters, which lead to the regimes described above. 
Such analysis must be based on the solution of the char- 
acteristic equation for Lq, which for N — 3 with trio = 
is: 



(30a) 



Now we can use Cardano formulas to determine the roots 
■ For each concrete set of soliton parameters it is easy 
to evaluate Cfc and to determine the asymptotic regime 
predicted by CTC. The complete analytic analysis for 
generic complex valued p and q is rather lengthy. For 
brevity and simplicity we limit ourselves to two special 
choices of the initial amplitudes: 



A) d 2 = 0, di = -d 3 = 2Av 0l 



B) 



di 



2A2/0 



-di 



(32a) 



(32b) 



In case A) it is possible to adjust the set of initial 
phases of the solitons so that p and q will take real values: 



Phi 



Ph 2 



{4,o}Li = {0,^2,o,0}, (33a) 



{S kfi }l =1 = {0, £2,0,20)2,0}, (33b) 



Pha 



{fc,o}Li 



{CU 2 ,o,2<5 2 ,o-^}, (33c) 



For real q and p the properties of the roots are determined 
by the sign of the discriminant Q of (30a) fl34| : 



Q 



(34) 



Skipping the details, we find that each of the regimes 
occurs for the following choices of the soliton parameters. 
For regimes (i) and (iii) we will study only case A); regime 
(ii) will be studied exhaustively for both cases A) and B). 

Regime (i) for real-valued p and q requires 
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dj) 



l 2' 



1 = \ ( a i rf3 + a 2 d i) + Q^d 1 d 2 d 3 , 



(30b) 



(30c) 



where d k = 2(v kfl - v ), v G = YJk=o V kfl/N is the aver- 
age amplitude. For the initial sets of soliton parameters 
specified above we find: 



a*; = ie e 



i{&k,o — <5fc+i,o)/2 



e = v e 



(31a) 



(31b) 



and 



p = el (e-^.o + e i(&.o-fc.o)) + _ + d 2 + rf 2) (31c) 



l A (d 3 e- i5 *'° + d^ 8 ^- 8 ^) + ^d 1 d 2 d 3 . (31d) 



Q <0. 



(35a) 



Obviously this requires p < 0. The condition ( |35a| ) is 
satisfied for each one of the following sets of parameters: 

i.a) Ai/ = 0, Phi with § < S 2 < %S 

i.b) Phi with 5 2 — tt and 



|A^ | < Ai>, 



cr,3; 



(35b) 



where Az/ cr 3 is defined by (23c ). 

i.c) |Az/o| > 0, Ph 2 with 5 2t o ^ and w. Note that in 
this and in the next subcases q = but p is complex. 

i.d) A^ = 0, {4,o}Li = {0,<$2,o, ^3,0} with £3,0 ^ 
and S 2 . -^ ^ § and \. 

Regime (ii) - three solitons form a bound state if 



all the three roots of the characteristic polynomial (30a) 
are purely imaginary, i.e. = iKk- Then by changing 
C = *Ci we § e t a polynomial 



Ci -pCi+*« = o 



(36) 



with the purely real roots Kk', therefore its coefficients p 
and iq must also be real and the well known classification 
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of the solutions of cubic equations (see |3J]) solves the 
problem exhaustively. In terms of the coefficients of j36| ) 
this situation takes place when 



Qi 



-Q < 0, p>0. 



(37) 



In terms of the soliton parameters ( J37| ) is satisfied in the 
following cases. 



Case A) (32a). For this choice of initial amplitudes 
we get soliton bound states in the following subcases: 
ii.a) Ph= {0,0,0} and \Au \ > 0; 
ii.b) Ph= {0,±tt,0} and 



\Au \ > 2y/2e = Av, 



cr,3i 



(38a) 



ii.c) Ph = {0, 7To, ttq} and Ph = {0, 0, 7To} where ttq = 
±7r and 



\Ais \ > 2 3 3/4 



eo; 



(38b) 



Case B) ( |32b ). For this choice of initial amplitudes 
we get soliton bound states in the following subcases: 
ii.d) Ph= {0,0,0} and \Au \ > 0; 
ii.e) Ph = {0, ±tt, 0} and 



\Av \ > 4V2e ; 



(38c) 



ii.f) Ph = {0,£ 2 ,o,0} and \Av \ > if cos(5 2 ,o > 0; if 
cos ($2,0 < then the bound state regime holds for 



\Au \ > 4V2 v /-cos(S2,oeo; 



(38d) 



As one may expect, the asymptotic behavior depends 
very much on the initial choice of phases. Thus in the 
cases ii.a) and ii.d) the bound state regime is obtained 
for any value of Avq , while in all other cases this regime 
is entered only if \Ai/q \ is larger than some critical value, 
which of c ours e, depe n ds on the i nitia l parameters, com- 
pare e.g. (§8a|), (§8b|), (§8cj) and (§8d|). 

Regime (iii) requires one of the two conditions: 



Q>0, q^0, 



= 0, 



i.e., 



P<0, 



±2 



(39a) 



(39b) 



(39c) 



For the soliton parameters (39a), ( 39b ) respectively 
mean: 

iii. a) Phi with § < 8 2 ,o < and Av Q = 

Av CT ,3y/— cos 62ft- Then p = and q is real. 

iiiib) \Avo\ > 0, Phi with <5 2 ,o = f ■ 

iii.c) \Av \ > A^ ci! 3, Phi with <5 2 ,o 7^ 0, n. 

We also note several critical sets of soliton parameters 
when one regime switches into another one. For example 



for Av = and Phi with S 2 ,o = § (or <5 2 ,o = if ~ 
— 5) regime (i) switches over into regime (ii). The same 
happens also for \Av \ = Av> cl 3 and Phi with <5 2 ,o = 7r - 

Th e complete study of the characteristic equation 
( |30a| ) in the general case, when both p and q are complex 
and nonvanishing will be given elsewhere. 

These three main types of asymptotic behavior exist 
also for cases with N > 3 particles. Then it is possible 
to choose the parameters and in such a way that 
the corresponding solution of the CTC will describe N- 
particle bound state; if, in addition, all the ratios (77/. — 
Vj)/{Vk~Vm) are rational, the corresponding solution will 
be periodic in time. One may also find the conditions 
under which this solution will develop singularities for 
finite t (like Ti = or 7 3j0 = above). 

Obviously, there exist also all intermediary types of 
asymptotic. For example, if we have K\ < . . . < — 
Kfe+i = • • • = Kk+m—i < ■ ■ ■ < k-n then we will get in the 
asymptotic an m-particle bound state and all the other 
particles will have free-motion asymptotic. 



III. COMPARISON BETWEEN CTC AND NLS 

In this section we compare the analytical results ob- 
tained in the Sections II with the numerical solutions of 
the unperturbed NLS (R[u] = 0) (0) with initial condi- 
tions (0). 

We have collected in tables, variety of runs. We will 
consider separately cases with different N and different 
type of asymptotic regimes. The quantity showing how 
well the solution of the CTC fits the NLS solution is given 
by the root mean square: 



\ 



Ni 

E 

a=l 



NLS 



a(i Q ) CTC ) 2 



Ni 



(40) 



By N\ above we have denoted the number of experimen- 
tal points t a , at which the numerical values are calcu- 
lated. 

In the case of regime (iii) we have also compared the 
mean values of the asymptotic velocities of the solitons 
calculated from NLS and CTC. 

N = 2. Initial conditions assuring different type of 
regimes were described in Sec. H 

Results comparing the NLS and the CTC are collected 
in table I and the first rows of tables III and IV. It is 
seen that the agreement is very good and improves with 
the increasing of A^o and r$. Transition from regime 
(i) to regime (ii) is also shown in table I. In addition to 
these results we show that the factor A\ (|l7|), see table 
IV, has very low values, which is a reason to call such 
regime quasi-equidistant. 

Degenerated case. It is realized when the two eigen- 
values of Lq are equal. This happens when Avq — Ai/ cr . 2 



and #2,0 = t- This case is very difficult for observing nu- 
merically, however our data show that the value of Av„^ 
is quite well predicted by CTC. 

Remark. We mention again that under initial condi- 
tions Sq — and Av = the CTC possesses singular 
solutions. It is known that the iV-soliton solution of (|l|) 
is always analytic, i.e. does not possess any singularities. 
Our numerical checks show that such sets of initial con- 
ditions correspond to solitons colliding (for small ro < 8) 
or coming very close to each other (for ro > 8). In the 
regions, where this happens, the adiabatic approxima- 
tion is no longer valid and the CTC does not match with 
the numeric solutions. For N — 2 our numerical cal- 
culations show that the first coalescence takes place at 
t = T2 + |ai/ =o/2 and the next ones tend to repeat peri- 
odically with period very close to T 2 + \au =o (HI)- Similar 
behavior is also found for N = 3 with T£~\a„ 0= o replaced 
by T^~\au =o- For N > 3 the picture quickly becomes 
rather complicated. 

Our conclusion is that for such initial conditions the 
quasi-equidistant propagation of the solitons is main- 
tained for a restricted period of t < T' 2 + |ai/ =o/2. 

N = 3. Numerical data for regime (i) for Az^o = 
and So close to n was given in |,|,||,|6|. In Table II we 
have collected data from runs for a variety of values for 
#2,o- From it we see a very good agreement between CTC 
and NLS when $2,0 is far from -|; this value is critical for 
a transition from regime (i) to regime (ii). In the vicinity 
of J the agreement is not too good, see case 3|7|0|^. It 
is also seen that agreement is better for ro = 8 than for 
ro = 7. 

Regime (ii) — bound state solitons. Data collected 
in table III shows very good agreement between NLS and 
CTC. The values of Ak collected in table IV for are rather 
low, especially for 2A^o = 0.15 and 0.20 and again we 
may call such propagation quasi-equidistant. 

As in the N = 2 case, the increasing of A^ and ro 
leads to better agreement between NLS and CTC (see 
Table III), and also leads to lower value of Ak (see Ta- 
ble IV). Another fact illustrated in these tables is the 
better agreement for So = than for So = n. 

Regime (iii). This type of regime is only available 
when N > 2. In this case two of the solitons form 
(nearly) bound state and the third one has free prop- 
agation. Examples of initial conditions leading to such 
type of regime are iii. a), iii.b) and iii.c). In table V 
we show the quantity 8\k as a criterion of agreement be- 
tween NLS and CTC. The soliton parameters for these 
runs correspond to the case iii.c) and regime (iii) must 
take place for all values of #2,0 7^ and tt. 

Another possibility to compare CTC and NLS is given 
in table VI where we show the mean values of the asymp- 
totic velocities of each soliton. We see that CTC always 
shows that two of the solitons form a bound state. This 
is not always what NLS shows, however the asymptotic 
values of the velocities of 2-nd and 3-rd solitons are quite 
close, so that we could consider it as a "quasi-bound" 
state. 



If Avo = A^ cr .3 and #2,0 — tt then all eigenvalues of Lq 
are equal and thus a degenerated case is realized. It is 
again very difficult to observe it numerically. The CTC 
gives logarithmically growing distance between solitons 
i.e. unbound state. 

N > 3. The cases listed above can also be extended 
to larger number of solitons showing the same type of 
regimes; of course now there are more mixed regimes pos- 
sible. To analyze them one should solve the characteristic 
polynomial of Lq and determine the soliton parameters, 
for which two or more of the «j's are different (regime (i)) 
or equal (regimes (iii) and (ii)). More details about them 
will be published elsewhere. Cases when the spectrum 
of L becomes degenerate, i.e. when two (or more) of the 
eigenvalues £jfc become equal, deserve special attention 
and further investigations. 

Here we concentrate only on regime (ii) as the one 
which might be important in long range fiber optics com- 
munications. In tables III and IV we compare CTC and 
NLS for soliton trains iV|8|2A^o|fc with k = and 1, and 
N = 2 to 5. The results listed in the table IV are from 
run lengths equal to t — 300 in dimcnsionless units; for 
regime (ii) we find periodic behavior with periods much 
less than 300 and small amplitudes of oscillations given 
by Ak. On several occasions, we extended the runs to 
lengths of t = 600 and 996, the same type of behav- 
ior and the same results for Ak- Thus we conclude that 
these sets of initial parameters lead to a very stable be- 
havior. Indeed, taking the amplitudes of the solitons to 
be increasing by the same amount of A^o we find for 
2A^o = 0.2 and ro = 8, an equidistant behavior with 
error A k < 10%. 

Regime (ii) is characteristic also for such trains with 
N > 5. However in such cases where one has to consider 
soliton trains with large number of solitons, a configu- 
ration with continually increasing amplitudes becomes 
impractical. Therefore we explore other possibilities for 
sets of initial soliton amplitudes. 

In fact, the first step towards the solution of this prob- 
lem was to use a soliton train of in-phase solitons with 
alternating amplitudes 2vi = 2v% = 1.0 and 2^2 = 2^4 = 
1.25 in Q and verified experimentally by ||. In Q the 
run length was equal to 141. This idea was experimen- 
tally verified in ||| where 20 Gbit /s single channel soliton 
transmission over 11 500 km using alternating amplitude 
solitons was reported. 

To check this idea we made a series of runs of N 
in-phase soliton trains with alternating amplitudes and 
N = 3 to 8, see Table VII. What we find is that the 
quasi-equidistant propagation takes place for t < T qc d 
with T^f ~ 250. For larger values of t, some of the soli- 
tons come rather close to each other. CTC also predicts 
the same type of behavior but with a larger value for 
T^ c ~ 460. This is the reason why, in Table VII, the 
numerical data from NLS shows non equidistant prop- 
agation for a run length of 300. We extended some of 
these runs to lengths of 660 and 996. The results for 
N = 3 show a structure close to a periodic one with a 
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period of about 2T^ S . So if the goal is to achieve a 
quasi-equidistant propagation to lengths t < then 
such configurations can be used. As an illustration of 
this fact, we have provided Fig. 2 where a train of 8 soli- 
tons with alternating amplitudes is shown. We see, that 
in the regions where the solitons tend to coalesce, the 
match between CTC and NLS is on a qualitative level. 

As one can see from Figs. 2-5, there are two distinct 
time scales in these figures. First, there are the oscil- 
lations with periods on the order of 10 or so. But at 
the same time, one notes that there is in general a much 
longer time scale, one that is on the order of 250 or larger. 
It is the motions on this longer time scale that violate the 
equidistant propagation, see Fig. 2. Note that in Fig. 4, 
when one can set the initial values so that the very long 
periodic motion is absent, then one has a very stable and 
equidistant motion. For larger run lengths, other config- 
urations must be used, see also 0. 

The next possibility which we explored was a "saw- 
like" configuration for the soliton amplitudes with three 
and four different amplitudes. The results are collected 
in Table VIII for N — 4 to 8. We see that configurations 
with only three different amplitudes show basically the 
same behavior as the alternating amplitudes case; the 
value of here is slightly larger, and again in the 

regions where the solitons tend to coalesce the match 
between CTC and NLS is only qualitative, see Fig. 3. 

Finally, for " saw-like" configurations of the amplitudes 
with four different amplitudes we get a substantially dif- 
ferent picture. It seems that here the value of is 
much larger than in the previous two cases. The CTC 
model shows a quasi-equidistant behavior with the same 
small values of Ak like in the last column of Table VIII 
to run lengths of 12 000. The numeric solution of NLS 
shows an excellent match with CTC to lengths of 1 200, 
see Fig. 4. 

IV. CONCLUSIONS 

A method for the description of the asymptotic behav- 
ior of the A^-soliton pulse trains of the NLS equation is 
proposed, based on the CTC model for the soliton inter- 
action. It describes correctly several qualitatively differ- 
ent classes of asymptotic regimes. Several sets of soliton 
parameters have been described for which the propaga- 
tion is quasi-equidistant. Such behavior with a conve- 
niently low value for Ak can be achieved by: a) taking 
soliton trains with At'o large enough. Increasing A^ 
too much may lead to violation of the condition ([h")) 
of the adiabatic approximation. Somewhat surprisingly, 
this nevertheless leads to a better agreement between the 
CTC and the NLS; b) increasing the value 2v$ of the aver- 
age amplitude; and c) increasing the distance ro between 
the neighboring solitons. 

The critical values of the soliton parameters, for which 
one regime switches over to another one have been eval- 



uated. We find that near these critical values, the match 
between CTC and NLS becomes worse, as one would ex- 
pect. 

These results have a natural generalization also for 
N > 3 soliton trains. The CTC-model provides one 
with a tool for constructing sets of initial data for the 
V-soliton trains which will possess a given asymptotic 
trait, which is also determined by the eigenvalues of 
L . 

The monotonically increasing amplitudes cannot be 
used in case one wants to send trains with larger num- 
ber of solitons. That is why we investigated also other 
possibilities for amplitudes for which the propagation is 
quasi-equidistant: alternating and "saw-like" configura- 
tions with three and four different values for 2vt-. We 
find that the alternating amplitudes configuration and 
the "saw- like" one with three different amplitudes may 
provide quasi-equidistant propagation up to run lengths 
on the order of about 230 or 300 correspondingly. The 
"saw-like" configurations with four different amplitudes 
show quasi-equidistant behavior to even much larger run 
lengths. 
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TABLE I. Comparison between NLS and CTC for each 
soliton (8xk)'- transition from regime (i) to regime (ii). Run 
length equals 300. 







Regime (i) 




Regime 


(ii) 


k 


2|7|08|1 


2|7|10|1 


2|7|12|1 


2|7|14|1 


2|7|16|1 


1 


0.31 


0.56 


1.58 


0.38 


0.22 


2 


0.54 


0.84 


1.90 


0.58 


0.36 




2|8|4|1 


2|8|6|1 




2|8|8|1 


2|8|10|1 


1 


0.040 


0.11 




0.35 


0.093 


2 


0.16 


0.29 




0.44 


0.16 



TABLE II. Comparison between NLS and CTC for each 
soliton (8xk)- regime (i). Run length equals 300. 



k 


3|8|0|li 


3|8|0|i| 


3|8|0|il 


3|8|0|1 


1 


0.89 


0.23 


0.16 


0.10 


2 


0.00 


0.00 


0.00 


0.00 


3 


0.89 


0.23 


0.16 


0.10 




3|7|0|M 


3|8|0|& 


3|8j0|£ 




1 


2.29 


0.68 


0.42 




2 


0.00 


0.00 


0.00 




3 


2.29 


0.68 


0.42 
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TABLE III. Comparison of NLS and CTC for each soliton 
(5xk)'- regime (ii).Run length equals 300. 



k 


2|8|10|0 


2|8|10|1 


2|8|15|0 


2|8|15|1 


2|8|20|0 


2|8|20|1 


1 


0.054 


0.093 


0.035 


0.040 


0.024 


0.027 


2 


0.061 


0.16 


0.050 


0.079 


0.044 


0.059 




3|8|10|0 


3|8|15|0 


3|8|15|1 


3|8|20|0 


3|8|20|1 




1 


0.098 


0.038 


0.084 


0.011 


0.041 




2 


0.21 


0.15 


0.25 


0.12 


0.19 




3 


0.12 


0.041 


0.066 


0.005 


0.035 






4|8|10|0 


4|8|15|0 


4|8|15|1 


4|8|20|0 


4|8|20|1 




1 


0.11 


0.031 


0.22 


0.021 


0.15 




2 


0.29 


0.21 


0.53 


0.22 


0.46 




3 


0.20 


0.12 


0.14 


0.079 


0.090 




4 


0.070 


0.041 


0.061 


0.029 


0.032 






5|8|10|0 


5|8jl5|0 


5|8|20|0 








1 


0.13 


0.048 


0.043 








2 


0.39 


0.30 


0.35 








3 


0.26 


0.19 


0.15 








4 


0.11 


0.061 


0.044 








5 


0.063 


0.055 


0.033 









TABLE IV. Comparison between NLS and CTC for each 
pair of neighbor solitons: value of Ak, ro = 8. Run length 
equals 300. 



k 


NLS CTC 


NLS CTC 


NLS CTC 




2|8|10|0 


2|8|15|0 


2|8|20|0 


1 


0.055 0.054 


0.029 0.026 


0.019 0.016 




3|S| 10 


3|8|15|0 


3|8|20|0 


1 


0.09 0.03 


0.049 0.014 


0.031 0.008 


2 


0.099 0.03 


0.053 0.014 


0.020 0.008 




4|8|10|0 


4|8|15|0 


4|8|20|0 


1 


0.12 0.09 


0.075 0.039 


0.058 0.019 


2 


0.14 0.13 


0.085 0.051 


0.064 0.026 


3 


0.08 0.09 


0.034 0.039 


0.020 0.019 




5|8|10|0 


5|8|15|0 


5|8j20|0 


1 


0.16 0.09 


0.11 0.04 


0.09 0.019 


2 


0.20 0.063 


0.13 0.026 


0.11 0.013 


3 


0.11 0.063 


0.059 0.026 


0.04 0.013 


4 


0.055 0.09 


0.02 0.04 


0.009 0.019 



TABLE V. Comparison between NLS and CTC for each 
soliton (5xk)'- regime (iii). Here amplitudes have the 
values: m = (0.9146,1.0,1.0854) (Av = Ai/ cr , 3 |r =7), 
n 2 = (0.9482,1.0,1.0518) (Av = Av CI , 3 \ ro =&)- Run length 
equals 300. 



k 


3|7|m|i 


3|7|m|§ 


3|7|m|i 


3|7|m|| 


3|7|m|| 


1 


0.24 


0.26 


0.68 


1.03 


1.16 


2 


0.83 


0.90 


0.68 


0.25 


0.84 


3 


0.28 


0.40 


0.50 


0.25 


0.13 




3|8|n 2 |i 


3|8|n 2 |i 


3|8|n 2 |i 


3|8|n 2 || 


3|8|n 2 |f 


1 


0.22 


0.12 


0.23 


0.35 


0.38 


2 


0.27 


0.28 


0.18 


0.23 


0.47 


3 


0.14 


0.19 


0.24 


0.13 


0.063 



TABLE VI. Comparison between NLS (mean asymptotic 
velocities) versus — 4£ fe from CTC and for each soliton: regime 
(iii) . Notations m , n 2 here mean the same as in the table 
above. Run length equals 300. 



k 


NLS 


CTC 


NLS 


CTC 


NLS 


CTC 




3|7|m 


4 


3|7|m 




3|7|m 




1 


-0.034 


-0.032 


-0.045 


-0.042 


-0.063 


-0.058 


2 


0.011 


0.016 


0.015 


0.021 


0.024 


0.029 


3 


0.018 


0.016 


0.024 


0.021 


0.032 


0.029 




3|7|m 


i 


3|7|m 


'A 

4 






1 


-0.073 


-0.067 


-0.074 


-0.067 






2 


0.032 


0.033 


0.037 


0.033 






3 


0.035 


0.033 


0.031 


0.033 








3|8jn 2 


4 


3|8|n 2 




3|8|n 2 




1 


-0.020 


-0.020 


-0.026 


-0.026 


-0.037 


-0.036 


2 


0.008 


0.010 


0.011 


0.013 


0.016 


0.018 


3 


0.011 


0.010 


0.014 


0.013 


0.019 


0.018 




3|8jn 2 


•A 


3|8|n 2 


'A 

4 






1 


-0.042 


-0.040 


-0.043 


-0.040 






2 


0.020 


0.020 


0.021 


0.020 






3 


0.021 


0.020 


0.020 


0.020 
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TABLE VII. Comparison between NLS and CTC for each 
pair of neighbor solitons: value of At for ro = 8. Amplitudes 
are ordered alternatingly; by "si" we have denoted trains with 



2ui 


= 1.0 and 2v2 = 


1.25 and "Is' 


means 


2i/i = 


1.25 and 




= 1.00. Run length equals 300. 








k 


NLS CTC 


NLS 


CTC 


NLS 


CTC 




3|8|sl|0 


3|8|ls|0 


4|8|sl|0 


1 


0.62 0.08 


0.012 


0.026 


0.33 


0.005 


2 


0.62 0.08 


0.012 


0.026 


0.42 


0.078 


3 








0.31 


0.005 




4|8|ls|0 


5|8|sl|0 


5 1 8 1 Is 1 


1 


0.31 0.005 


0.23 


0.014 


0.59 


0.006 


2 


0.42 0.078 


0.03 


0.058 


0.59 


0.031 


3 


0.33 0.005 


0.03 


0.058 


0.59 


0.031 


4 




0.23 


0.014 


0.59 


0.006 




6|8|sl|0 


6|8|ls|0 


7|8|sl|0 


1 


0.24 0.005 


0.28 


0.005 


0.26 


0.009 


2 


0.024 0.039 


0.27 


0.039 


0.14 


0.05 


3 


0.036 0.004 


0.036 


0.004 


0.10 


0.006 


4 


0.27 0.039 


0.024 


0.039 


0.10 


0.006 


5 


0.28 0.005 


0.24 


0.005 


0.14 


0.05 


6 








2.04 


0.073 




7|8|ls|0 


8|8|sl|0 


8 1 8 1 Is 1 


1 


0.27 0.005 


0.25 


0.005 


0.31 


0.005 


2 


0.27 0.035 


0.14 


0.04 


0.31 


0.04 


3 


0.026 0.002 


0.10 


0.003 


0.13 


0.003 


4 


0.026 0.002 


0.13 


0.005 


0.13 


0.005 


5 


0.27 0.035 


0.13 


0.003 


0.10 


0.003 


6 


0.27 0.005 


0.31 


0.04 


0.14 


0.04 


7 




0.31 


0.005 


0.25 


0.005 



TABLE VIII. Comparison between NLS and CTC for each 
pair of neighbor solitons: value of At . The train is a saw- like 
one with amplitudes 3sw — ► {0.85,1.00,1.15,0.85,...} and 
4sw {0.85,1.00,1.15,1.30,0.85,...}. Run length equals 
1200. 



k 


NLS 


CTC 


NLS 


CTC 


NLS CTC 




4|8|3sw|0 


5|8|3sw|0 


5|8|4sw|0 


1 


0.19 


0.07 


0.47 


0.07 


0.05 0.04 


2 


0.09 


0.05 


0.25 


0.06 


0.06 0.04 


3 


1.13 


0.04 


0.39 


0.05 


0.02 0.04 


4 






0.41 


0.06 


0.01 0.01 




6|8|3sw|0 


6|8|4sw|0 


7|8|3sw|0 


1 


0.46 


0.20 


0.10 


0.03 


0.79 0.30 


2 


0.23 


0.22 


0.06 


0.04 


0.61 0.44 


3 


0.36 


0.12 


0.02 


0.04 


0.47 0.22 


4 


0.40 


0.22 


0.05 


0.02 


0.66 0.35 


5 


0.15 


0.20 


0.08 


0.02 


0.54 0.45 


6 










0.18 0.18 




7|8|4sw|0 


8|8|3sw|0 


8|8|4sw|0 


1 


0.09 


0.03 


0.84 


0.19 


0.09 0.03 


2 


0.06 


0.04 


0.38 


0.30 


0.06 0.03 


3 


0.02 


0.03 


0.26 


0.13 


0.02 0.03 


4 


0.05 


0.02 


0.36 


0.28 


0.05 0.02 


5 


0.08 


0.03 


0.22 


0.33 


0.08 0.03 


6 


0.05 


0.03 


0.49 


0.06 


0.05 0.03 


7 






0.52 


0.15 


0.02 0.03 
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FIG. 1. Plot of the functions Ai-j, m (Auo) in Eq. 
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FIG. 2. Propagation of 8 solitons with alternating amplitudes 2^ 2 fc-i = 1.0, 2v 2 k = 1.25, r = 8. 
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FIG. 3. Propagation of 8 solitons with a "saw"-like configuration of the amplitudes: 2i/ 3fc _2 — 0.85, 2^ 3fe _i = 1.0, 2v; ik — 1.15, 
r = 8. 
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FIG. 4. Propagation of 8 solitons with a "saw"-like configuration of the amplitudes: 2v4ks = 0.85, 2i/4k-2 = 1-0, 
2^4fc-i = 1.15, 2vt,k = 1-30, r = 8. 
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FIG. 5. Propagation of 3 solitons with alternating amplitudes: 2^1,0 = 2^ 3i o = 1.00, 2f2,o = 1.25, ro = 8. 
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